Can the nearly conformal sextet gauge model hide the Higgs impostor? 
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Abstract 

{N) New results are reported from large scale lattice simulations of a frequently discussed strongly interacting gauge theory with 
t-H a fermion flavor doublet in the two-index symmetric (sextet) representation of the SU(3) color gauge group. We find that the 
< ^ chiral condensate and the mass spectrum of the sextet model are consistent with chiral symmetry breaking in the limit of vanishing 
fermion mass. In contrast, sextet fermion mass deformations of spectral properties are not consistent with leading conformal scaling 
Mh behavior near the critical surface of a conformal theory. A recent paper could not resolve the conformal fixed point of the gauge 
^ coupling from the slowly walking scenario of a very small nearly vanishing yS-function Q. It is argued that overall consistency with 
our new results is resolved if the sextet model is close to the conformal window, staying outside with a very small non-vanishing 
jS-function. The model would exhibit then the simplest composite Higgs mechanism leaving open the possibility of a light scalar 
| state with quantum numbers of the Higgs impostor. It would emerge as the pseudo-Goldstone dilaton state from spontaneous 
symmetry breaking of scale invariance. We will argue that even without association with the dilaton, the scalar Higgs-like state can 



be light very close to the conformal window. A new Higgs project of sextet lattice simulations is outlined to resolve these important 



q, questions. 
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^ 1. Introduction 

Q\ The stunning discovery of the 125 GeV Higgs-like parti- 
es cle at the Large Hadron Collider [1 2 1 does not exclude new 
physics beyond the Standard Model (BSM) in the framework 
Q\ of some new strongly-interacting gauge theory with a compos- 
ite Higgs mechanism, an idea which was outside experimental 
reach when it was first introduced as an attractive BSM sce- 
. . nario B HTTI . The original framework has been considerably 
!>■ extended by new explorations of the multi-dimensional the- 
k> ory space in fermion flavor number, the choice of color gauge 
j_j group, and fermion representation lfT2l - l2"Tl . Systematic and 
non-perturbative lattice studies play an important role in stud- 
ies of this extended theory space [22 -62 1. Even without spin 
and parity information, the new Higgs-like particle with decay 
modes not far from that of the Standard Model brings new focus 
and clarity to the search for the proper theoretical framework. 

One example is the light dilaton as a pseudo-Goldstone par- 
ticle of spontaneous breaking of scale invariance that has been 
featured in recent phenomenological discussions as a viable in- 
terpretation of the discovery [63, 64 1. Nearly conformal gauge 
theories serve as theoretical laboratories for realistic implemen- 
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tations of this scenario [65 -73]. Unfortunately, a credible real- 
ization of the idea as a strongly interacting BSM gauge theory 
is still lacking. We investigate here a candidate theory with a 
fermion flavor doublet in the two-index symmetric (sextet) rep- 
resentation of the SU(3) color gauge group close to the confor- 
mal window, if it can hide a light Higgs-like scalar state with or 
without dilaton-like interpretation. 

The sextet force and a new fermion doublet driving elec- 
troweak symmetry breaking was introduced in QCD a long time 
ago by Marciano fl4l . Early pioneering lattice work, limited to 
the quenched approximation at that time, investigated the sextet 
fermion representation ifTBTl . The main difference in the model 
we investigate here is the introduction of a new SU(3) gauge 
force not associated with QCD gluons and motivated by ideas 
of compositeness from a new super-strong force. After chiral 
symmetry breaking we find three massless Goldstone pions in 
the spectrum providing the minimal realization of the Higgs 
mechanism, just like in the original technicolor idea (4] [5] . 
The important new ingredient is the sextet representation of the 
fermion doublet which brings the model very close to the con- 
formal window as indicated in a recent paper |3 1. The accuracy 
of the very small nearly vanishing /3-function in difficult sim- 
ulations could not resolve the existence of a conformal fixed 
point gauge coupling from the alternative slowly walking sce- 
nario. When combined with our observation of chiral symmetry 
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breaking OfSB) reported here for small fermion mass deforma- 
tions, the overall consistency of all simulations is resolved if 
the sextet model is close to the conformal window with a very 
small non-vanishing /J-function (see, also l50ll5ll ). In this case 
the model exhibits the simplest composite Higgs mechanism 
and leaves open the possibility of a light scalar state with quan- 
tum numbers of the Higgs impostor emerging as the pseudo- 
Goldstone dilaton state from spontaneous symmetry breaking 
of scale invariance. Even if scale symmetry breaking is entan- 
gled with^SB without dilaton interpretation, a light Higgs-like 
scalar state can emerge from the new force close to the confor- 
mal window. Our new Higgs project with lattice simulations in 
the sextet model may resolve these important problems. 

In section 2 we will outline the computational strategy and 
the simulation set-up including the important treatment of fi- 
nite size effects. In section 3 results on the chiral condensate are 
presented with extrapolation to the massless fermion limit. The 
spectrum is presented in Section 4 and compared with the ^ SB 
hypothesis. In Section 5 it is shown that fermion mass deforma- 
tions of spectral properties are not consistent in the model with 
conformal scaling behavior near the critical surface of a con- 
formal theory. Section 6 will describe the new Higgs project 
to determine the scalar ++ mass spectrum when disconnected 
diagrams are included in the calculations. Closely related to the 
dilaton interpretation, we also outline in Section 6 the important 
role of the non-perturbative gluon condensate and our strategy 
for investigating it within our new Higgs project. 

2. Computational strategy and lattice simulations 

Probing ;fSB, and conformal behavior for comparison, we 
extrapolate the spectrum to infinite volume at fixed fermion 
mass m. In large volumes the leading finite size corrections 
are exponentially small and dominated by the lowest state of 
the spectrum which has pion quantum numbers. From the mass 
spectrum, extrapolated to infinite volume, we can probe the pat- 
tern of ^ SB when small fermion mass deformations are simu- 
lated close to the massless limit. We also probe the hypothesis 
of mass deformed conformal scaling behavior. Our results, as 
we report here, strongly favor the^SB hypothesis. 

2.1. The algorithm and simulation results 

We have used the tree-level Symanzik-improved gauge ac- 
tion for all simulations reported in this paper. The conventional 
P = 6/g 2 lattice gauge coupling is defined as the overall factor 
in front of the well-known terms of the Symanzik lattice action. 
Its values are j3 = 3.20 and /3 = 3.25 for our simulations. The 
link variables in the staggered fermion matrix were exponen- 
tially smeared with two stout steps [74|; the precise definition 
of the staggered stout action was given in ll75ll . The RHMC al- 
gorithm was deployed in all runs. The fermion flavor doublet 
requires rooting in the algorithm. For molecular dynamics time 
evolution we applied multiple time scales l76l and the Omelyan 
integrator |77|. Our error analysis of hadron masses used corre- 
lated fitting with double jackknife procedure on the covariance 
matrices |78 1. The time histories of the fermion condensate, the 



plaquette, and correlators were used to monitor autocorrelation 
times in the simulations. 

We have new simulation results at [5 = 3.2 in the 
fermion mass range m = 0.003 - 0.010 on 24 3 x 48, 28 3 x 
56, and 32 3 x 64 lattices. Five fermion masses at 
m = 0.003,0.004,0.005,0.006,0.008 are used in most fits. A 
very large and expensive 48 3 x 96 run was added recently at 
m = 0.003 to control finite size effects. We also have new sim- 
ulation results at/3 = 3.25 in the mass range m = 0.004 - 0.008 
on 24 3 x 48, 28 3 x 56, and 32 3 x 64 lattices. 

2.2. Finite size effects 

Infinite-volume extrapolations of the lowest state in the spec- 
trum with pion quantum numbers, the related F„, and the con- 
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Figure 1: Finite volume dependence at the lowest fermion mass for fi = 3.2. 
The form of ~gi(A, rfj is a complicated infinite sum which contains Bessel func- 
tions and requires numerical evaluation [79 ]. Since we are not in the chiral log 
regime, the prefactor of the ~g\(A,r}) function was replaced by a fitted coeffi- 
cient. The leading term of the function ~j>\ (A, rj) is a special exponential Bessel 
function Ki{A) which dominates in the simulation range. 



densate (if/ij/} are shown in Figure [T| where ~g\(A, if) describes 
finite volume corrections from the exchange of the lightest pion 
state with A = M„L and lattice aspect ratio rj = T/L, similarly 
to what was introduced in ll80l . The fitting procedure approxi- 
mates the leading treatment of the pion which wraps around the 
finite volume, whether in chiral perturbation theory Ctpt), or in 
Luscher's non-perturbative finite size analysis ||8T|| . This equiv- 
alence relaxes the requirement on the fitted parameters cm,cf,c\ 
to agree with 1-loop^PT as long as the pion is the lightest state 
dominating the finite volume corrections. It should be noted 
that the form of the fitting function ~gi(A, rj) does not commit to 
the chirally broken phase. At fixed fermion mass m, the leading 
exponential term of the function is also the expected behavior 
in the conformal phase with mass deformation. The asymptotic 
exponential form simply originates from the lightest state wrap- 
ping around the volume once emitted from and re-absorbed by 
the composite state whose sensitivity to finite volume correc- 
tions is being investigated. The analysis is therefore applicable 
to both mass deformed phases with different symmetry proper- 
ties. 

The infinite-volume limits of M„, F n , and (i/nff) for m = 0.003 
at B = 3.2 were determined self-consistently from the fitting 
procedure. Similar fits were applied to other composite states. 
The value of M„ in the fit of the top plot in Figure [T] was de- 
termined from the highly non-linear fitting function and used 
as input in the other two fits. Based on the fits at m = 0.003, 
the results are within one percent of the infinite-volume limit at 
M„L = 5. In the fermion mass range m > 0.004 the condition 
M„L > 5 is reached at L — 32. Although it will require high 
precision runs to test, we expect less than one percent residual 
finite size effects in the 32 3 x 64 runs for m > 0.004. Based 
on these observations, we will interpret the results from the 
32 3 x 64 runs for m > 0.004 as infinite-volume behavior in 
mass deformed chiral and conformal analysis. 

3. The chiral condensate 

Our simulations show that the chiral condensate (ifnfr) is con- 
sistent with ^SB and remains non-vanishing in the massless 
fermion limit. It has the infinite-volume spectral representation, 

<^> = -2m- f dA-^-, (1) 
Jo m z + A 1 

which is UV-divergent when the cutoff A is taken to infinity. 
The divergences are isolated by writing the integral of the spec- 
tral representation in twice subtracted form [82 1, 

W> = -2m • f dA 

Jo m 1 + A 1 

~ 2m5 ' f ~!x P } A \j + ci(fl) ■ m + c 3 (a) ■ m 3 . (2) 
J^ A 4 m 2 + A 2 

The first integral in Eq. Q isolates the infrared part and re- 
covers the well-known relation (i/fifr) = -np(0) in the m — > 
limit 1 83 1. The linear fermion mass term c\{a) ■ m is a quadrati- 
cally divergent UV contribution as cC 2 ■ m with lattice cutoff a. 



There is also a very small third-order UV term CT,(a)-m 3 without 
power divergences which is hard to detect for small m and has 
not been tested within the accuracy of the simulations. 

IR finite contributions to the condensate from the chiral La- 
grangian are connected at the low energy scale /i with the first 
integral in Eq. Q. In the chiral expansion of the condensate 
there is an m-independent constant term which is proportional 
to BF 2 , a linear term proportional to B 2 ■ m, a quadratic term 
~ B 3 F~ 2 ■ m 2 , and higher order terms, in addition to logarithmic 
corrections generated from chiral loops. The expansion in the 
fermion mass is expressed in terms of low energy constants of 
chiral perturbation theory, like B and F |84|. 
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Figure 2: The chiral condensate and its reduced form with subtracted derivative 
(both have to converge to the same chiral limit) are shown in the top plot with 
linear fit to the condensate. The data without derivative subtraction cannot de- 
tect higher order fermion mass terms with significant accuracy. The fit to the 
reduced form with subtracted derivative is defined in the text and shown in the 
magnified lower plot. A linear term is not included in this fit since the sub- 
tracted derivative form approximately eliminates it. The value of do at m = is 
shown to be consistent with the direct determination of cq from the chiral limit 
of (if/ift). The consistency is very reassuring since the two results are derived 
from independent determinations. For m = 0.003 the data from infinite-volume 
extrapolation were used in the fit. As we explained earlier, at higher m values 
the largest volume 32 3 X 64 runs were used for the condensate and its derivative 
subtraction. 

We used two independent methods for the determination of 
the chiral condensate in the massless fermion limit. In the first 
method fits were made directly to (tfn/r) with constant and linear 
terms in the fitted function. Quadratic and third order terms 
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are hard to detect within the accuracy of the data. The result 
is shown in the top plot of Figure [2] When the quadratic term 
is added to the fit, the massless intercept cq = (tf/ip) m= o from 
the quadratic fit agrees with the one from the linear fit and the 



quadratic fit coefficient in C2 ■ m is zero within fitting error. 

For an independent determination, we also studied the sub- 
tracted chiral condensate operator defined with the help of the 
connected part Xconn of the chiral susceptibility x, 



[i 



dm. 



(\j/\f/) - to ■ Xc 
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The derivatives d/dm and d/dm v are taken at fixed gauge 
coupling B. The derivative d/dm v is defined in the partially 
quenched functional integral of (4rt//) pq with respect to the va- 
lence mass to v and the limit to v = m is taken after differentia- 
tion. The removal of the derivative term significantly reduces 
the dominant linear part of the (ifnfr) condensate without chang- 
ing the intercept in the m = limit. Once the derivative term is 
subtracted, the first non-perturbative IR contribution, quadratic 
in to, is better exposed. The two independent determinations 
give consistent non-vanishing fit results in the massless chiral 
limit as shown in the lower plot of Figure [2] 

The independent determinations of the non-vanishing con- 
densate in the chiral limit with separate fits cq - (ijrift) m =Q and 
do = (ifn/r)m=o are consistent with each other but differ from the 
GMOR ED relation (ij/ip) = 2BF 2 by a factor of two. As shown 
in the next section, the value of 2B is determined in lattice units 
from the pion spectrum using the leading M 2 = 2B ■ to rela- 
tion. We find the numerical value 2B = 6.35(21) as shown in 
the top plot of Figure [3] F is determined from the pseudoscalar 
correlator which satisfies the PCAC relation. We find in lattice 
units the numerical value F = 0.0279(4) from the lower plot 
of Figure |3] with 2BF 2 = 0.0049(2). Both sides of the GMOR 
relation are sensitive to cutoff effects in B and F at bare lattice 
coupling B = 3.2. Our preliminary fits based on staggered chi- 
ral perturbation theory indicate that cutoff effects modifying the 
continuum values of B and F are likely sources of the discrep- 
ancy [86]. Some increase in the cutoff dependent values of B 
and F, which is the observed trend, would bring the two sides 
of the GMOR relation in agreement. 



test the two contrasting hypotheses, one with ;f SB and the other 
with conformal behavior. As a first step to address the removal 
of finite lattice spacing, we will compare the Goldstone and 
non-Goldstone pion spectra at two different lattice spacings to 
probe the restoration of taste symmetry for staggered fermions 
as the lattice spacing is decreased. 

4.2. The Goldstone pion and F„ 

The chiral Lagrangian describes the low energy theory of 
Goldstone pions and non-Goldstone pions in the staggered lat- 
tice fermion formulation. It will be used as an effective tool 
probing the ^SB hypothesis at finite fermion masses including 
extrapolation to the massless chiral limit. Close to the chiral 
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4. Spectral tests of the^SB hypotheses 

4.1. Strategy and challenges of the spectrum analysis 

Spectrum calculations in a gauge theory with massless 
fermions require important and difficult lattice extrapolations: 

(1) Extrapolation from finite lattice size to infinite volume, 

(2) Extrapolation to the massless fermion limit, 

(3) Extrapolation in lattice spacing to the continuum. 

All three issues will be addressed as we present details of the 
spectrum analysis in this section. The strategy of finite size 
corrections was explained in Section 2 and it will be applied 
here. Extrapolation from finite fermion masses will be used to 



Figure 3: Polynomial fits from the analytic mass dependence of the chiral La- 
grangian without logarithmic loop corrections are shown for the Goldstone pion 
and F„. The dashed line in the top plot for the Goldstone pion shows the leading 
linear contribution. 

limit, the pion spectrum and the pion decay constant F„ are 
organized in powers of the fermion mass to which is an input 
parameter in the simulations. Chiral log corrections to the poly- 
nomial terms are generated from pion loops [ 87 ] . Their analysis 
will require an extended dateset with high statistics. 

In Section 2 we presented results of infinite-volume extrap- 
olations. The effects are largest at to = 0.003 in our dataset 
and the infinite-volume limits of M n and F„ were shown for 
to = 0.003 for fixed lattice cutoff and bare coupling B = 3.2. 
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Similar fits were applied to the chiral condensate and compos- 
ite states in the spectrum at m = 0.003. Based on the analysis 
at m = 0.003, we determined that the infinite-volume limit is 
reached at M n L = 5 within one percent accuracy. It is expected 
that similar or better accuracy is reached for M„L > 5 at higher 
m values in all states of the spectrum. In the fermion mass range 
m > 0.004 the condition M„L > 5 is reached at L = 32. Based 
on these observations, in fits to the observed pion spectrum and 
F„ we will use infinite-volume extrapolation at m - 0.003 and 
treat the 32 3 x 64 runs for m > 0.004 as if the volume were 
infinite. 

In Figure [3] we used the local pion correlator with noisy 
sources to extract M„ and F„. The correlator is tagged as the 
PCAC channel since the PCAC relation, based on axial Ward 
identities, holds for this correlator and F„ can be directly deter- 
mined from the residue of the pion pole. The other staggered 
meson states and correlators we use are defined in [88]. For 
example, what we call the non-Goldstone scPion and the /o 
meson are identified in correlator I of Table 1 in [88) • Simi- 
larly, the non-Goldstone i5Pion is from correlator VII, the non- 
Goldstone ijPion is from correlator VIII, and the rho and Al 
mesons are from correlator III of Table 1 in ll88l . We measure 
the Goldstone pion in two different ways, with one of them de- 
fined above and the other is correlator II of Table 1 in |88|. 
For baryon states in the sextet fermion representation, not pre- 
sented here, we use our own construction of correlators which 
are different from the baryon correlators of 1 88 1 . 

Based on the analytic fermion mass dependence of the chi- 
ral Lagrangian, and using the lowest four fermion masses, good 
polynomial fits were obtained without logarithmic loop correc- 
tions as shown in Figure[3]for M„ and F n . Although we could fit 
M„ and F„ with the continuum chiral logarithms included, the 
two sets of F and B values from separate fits to M„ and F„ are 
not quite self-consistent. Rooted and partially quenched stag- 
gered perturbation theory is a useful procedure at finite lattice 
spacing for simultaneous fits of M„ and F„ with a consistent 
pair of F and B values |89 90 1. The explicit cutoff dependent 
corrections to the F and B parameters would require further 
testing at weaker gauge couplings and a set of valence fermion 
masses. 

We made the first step in this direction by adding a new run 
set to our database at B — 3.25. In Figure [4] we show taste- 
breaking effects in two pion spectra for comparison. We find 
significant reduction in taste breaking at smaller lattice spac- 
ing at the weaker coupling. Our staggered perturbation theory 
analysis will be presented in a longer follow-up report which 
will also include other results from the new runs at the weaker 
coupling /? = 3.25 H61. 

4.3. Taste breaking in the non-Goldstone pion spectrum 

The non-Goldstone pion spectra, quite different from the one 
found in QCD, are shown at B = 3.2 in the top plot of Fig- 
ure |4] using standard notation, introduced earlier. The non- 
Goldstone i5Pion is split from the Goldstone pion and remains 
exactly degenerate with the non-Goldstone scPion, a similar 
feature in QCD. The new feature is the mass dependence of 
the split between the Goldstone pion and the non-Goldstone 
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Figure 4: The top plot in the figure is the spectrum at /? = 3.2. It shows the 
polynomial fit of the Goldstone pion (magenta points). The red points are the 
non-Goldstone scPion data covering the green i5Pion data with complete de- 
generacy. The slightly split ijPion is shown with cyan color. The lower plot 
in the figure is the spectrum at fi = 3.25. In identical notation it displays the 
improvement in taste splitting with a considerably less taste-broken spectrum 
when plotted on the same scale. 



i5Pion with non-parallel slopes of the fitting functions. The 
non-Goldstone ijPion is further split from the i5Pion with a 
small mass-independent offset. Although taste breaking effects 
appear substantial on the scale of the plot, they are comparable 
with those from the HISQ action when the lattice spacings are 
matched ll9"T1 . The trends of the splits, particularly the fan-out 
structure and the lack of parallel equi-spaced splits with a con- 
stant slope determined by B is characteristic of gauge models 
as they get close to the conformal window. A very small resid- 
ual mass at m = is consistent with fits for the non-Goldstone 
pion states and decreases as we lower the lattice spacing with 
the weaker coupling at B = 3.25. This is shown in the lower 
plot of Figure [4] which exhibits a similar structure for the same 
pion states as the top figure but on a significantly more col- 
lapsed scale. Taste breaking is reduced considerably. It will be 
interesting to conduct a full analysis of all data on the finer lat- 
tice scale, closer to the continuum limit, and compare with the 
results presented here on the coarser lattice scale [86|. 
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4.4. The p and A i parity partner states 

It is useful and important to investigate the chiral limit of 
composite hadron states separated by a gap from the Goldstone 
and non-Goldstone pion spectra. The baryon mass gap in the 
chiral limit can provide further evidence for ^SB but our pre- 
liminary results are not shown here. Hadron masses of par- 
ity partners also provide important information with split parity 
masses in the chiral limit. This is particularly helpful not only to 
confirm^SB but to obtain a first estimate on the S parameter for 
probing the model against electroweak precision tests ||92ll . As 
an example, we will briefly review our results for the p meson 
state and its parity partner, the A; meson. Particularly interest- 
ing is the p — Aj mass splitting with parity violation. Figure [5] 
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Figure 5: Linear fit to the p meson mass is shown in the top plot of the figure. 
The lower plot shows the linear fit to the Ai meson superimposed on the p 
meson plot. The parity split is quite visible with varying size errors in the fitted 
m range. 

shows fits to the p meson and its A\ parity partner. The top 
plot is a linear fit to the p meson with a non-vanishing mass at 
m - 0, consistent with ^SB. The lower plot shows the linear 
fit to the A] meson. Both states extrapolate to non-vanishing 
masses in the chiral limit. The split appears to be significant 
for all fermion masses but the error is too large to resolve the 
chiral limit. More work with higher statistics is needed on this 
correlator before conclusive results can be obtained. 

5. Spectral tests of the conformal scaling hypothesis 

Under the conformal scaling hypothesis, the mass M n and the 
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Figure 6: The two plots represent separate conformal fits to M n (top) and F„ 
(bottom). The separate fits have reasonable^ 2 values but the incompatibility of 
the fitted y values disfavors the conformal hypothesis in its leading form. 

decay constant F„ are given at leading order by M n = c„-m t 1+y 
and F„ = Cf-m l ^ 1+y . The coefficients c„ and cp are channel spe- 
cific but the exponent y is universal in all channels B3TI461 . The 
leading scaling form sets in for small m values, close to the crit- 
ical surface. According to the hypothesis, there is an infrared 
conformal fixed point on the critical surface which controls the 
conformal scaling properties of small mass deformations. All 
masses of the spectrum can be subjected to similar conformal 
scaling tests, but we will mostly focus on accurate data in the 
M n and F„ channels. 

When M„ and F„ are fitted separately in the range of the four 
lowest fermion masses closest to the critical surface, we get 
reasonable x 2 values for the fits, as shown in Figure [6] How- 
ever, the incompatibility of the fitted y values disfavors the hy- 
pothesis, inconsistent with mass deformed conformal behavior. 
The conflicting simultaneous fits to universal conformal form 
with the same y for the Goldstone pion and the F„ decay con- 
stant are illustrated in Figure [7] Fitting to the pion mass sepa- 
rately requires y — 1 .040(73) while the separate F„ fit is forc- 
ing y = 2.20(15). In the combined fit they compromise with 
y = 1.53(28) and the unacceptable ^ 2 /dof of 44.5. It is impor- 
tant to note that the exponent y for the fit to M n only is what 
^SB would prefer. The separate conformal exponent y for F n 
is large to force to the origin the linear string of data which ex- 
trapolate to a finite constant in^SB. This creates conflict with 
the universal exponent y in the conformal analysis. 

From the tests we were able to perform, the sextet model is 
consistent with ^SB and inconsistent with conformal symme- 



6 



fit from combined conformal fit of and in PCAC channel 

M = C m 1 ' 1 "' 1 conformal fit p=3.2 
c = 1.41 ±0.33 




0.01 
0.008 
0.006 
0.004 
0.002 


j -0.002 
-0.004 
-0.006 
-0.008 
-0.01 



0.01 
0.008 
0.006 
0.004 
0.002 


j -0.002 
-0.004 
-0.006 
-0.008 
-0.01 



input from volumes 32 3 x 64, 48 3 : 
m fit range: 0.003 - 0.006 



0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.0 

m 

M fit residuals in combined fit 



fitting function: m 1 ' 






0=1.41 ±0.33 






V= 1.53 ± 0.28 




i 


X 2 /dof= 44.45 






- x 2 global sum= 222.2 


i 




X 2 sum (from M )=122 






I 








m fit range: 0.003 


- 0.006 



F fit from combined conformai fit of M and F in PCAC chann 




input from volumes 32 3 x 64, 48 3 x 96 
rtl fit range: 0.003 - 0.006 



F fit residuals in combined fit 



fitting function: c p m 1 ' 1 '*"' 




c p = 0.362 ±0.088 




y= 1.53 ±0.28 




X 2 /dof= 44.45 






* 


X 2 sum (from F_)= 100.1 




m fit range: 0.003 - 0.006 





Figure 7: The first plot shows the simultaneous conformal fit result for the pion 
mass, while the second displays the M„i residuals. The last two plots show the 
simultaneous fit result for the pion decay constant and the F„i residuals." The 
combined fit forces y = 1.53(28) with an unacceptable ^ 2 /dof of 44.5. 



try. It will require further investigations to show that subleading 
effects cannot alter this conclusion. We will consider compre- 
hensive conformal finite size scaling (FSS) tests which do not 
rely on infinite-volume extrapolation in the scaling fits. Confor- 



mal FSS was extensively applied to a different much discussed 
model with twelve fermion flavors in the fundamental repre- 
sentation of the SU(3) color gauge group E5ll . These kinds of 
tests are at a preliminary stage in the sextet project requiring 
new runs and systematic analysis. The FSS analysis of the ex- 
isting dataset of this paper, when smaller volumes are included, 
disfavors the conformal hypothesis similarly to what we just 
presented in the infinite-volume limit. It is difficult to reconcile 
^SB and large exponents in the fermion mass dependence with 
the low value of y defined by the chiral condensate using the 
Schodinger functional for massless fermions Q. 

6. The new sextet Higgs project 

If ^SB of the sextet model is confirmed in the massless 
fermion limit, its potential relevance for the realization of the 
composite Higgs mechanism is self-evident. The three Gold- 
stone pions of the model have the perfect match for providing 
the longitudinal components of the W ± and Z bosons. The re- 
maining most important issues are: (1) to calculate the mass 
of the ++ state when the disconnected part of correlator I in 
Table 1 of [88 1 is included; (2) the determination of the non- 
perturbative gluon condensate on the lattice to clarify the dila- 
ton connection if the Higgs particle turns out to be light; (3) a 
more precise determination of the running coupling for which 
we will deploy our new method based on the gradient flow of 
the gauge field in finite volume ||93l . We will outline in some 
details the first and second issues. 

6.1. The /o state in the ++ channel 

Figure[8]shows the fermion mass dependence of the fo meson 
without including the disconnected part of correlator I in Table 
1 of [88 1 . The non-Goldstone scPion and fo are parity partner 
states in this correlator. The quantum numbers of the fo meson 
match that of the ++ state in the staggered correlator. Close 
to the conformal window the fo meson is not expected to be 
similar to the <x particle of QCD. The full fo state including the 
disconnected diagrams could replace the role of the elementary 
Higgs and act as the Higgs impostor if it turns out to be light. It 
is very difficult to do the full calculation including the discon- 
nected diagram which is the main part of our next generation 
sextet Higgs project. First, we will discuss preliminary results 
which ignore the disconnected part. The challenges will be out- 
lined in the effort to include the disconnected part. 

The linear fit from the connected diagram is shown in Fig- 
ure [8] It has a non-zero intercept in the chiral limit with a mass 
more than five times F so it corresponds to a heavy state and 
not a Higgs candidate. Since the fo state is the parity part- 
ner of the non-Goldstone scPion in the full correlator, the two 
states would become degenerate in the chiral limit with unbro- 
ken symmetry. Close to the conformal window it is reasonable 
to expect that the disconnected diagram will dramatically re- 
duce the fo mass and its split from the scPion when the chiral 
limit is taken. This will leave the full fo state a viable Higgs 
candidate before new simulations resolve the issue and perhaps 
eliminate this attractive scenario. 
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effective mass M™ 11 from ++ connected correlator 

M to° n = M o + c i m l 1=3 " 2 
6 - meson mass requires missing disconnected part 
M = 0.1555 ±0.0070 

S °- 5 " C,= 25.3 ± 1.3 
8 ° 1 




m fit range: 0.003 -0.010 

I 1 1 1 

0.005 0.01 0.015 

m 

Figure 8: The linear fit is shown to the mass of the ++ fo meson from the 
connected part of correlator I in Table 1 of (.881. For comparison, the scPion 
which is the parity partner of the fo meson in the correlator is replotted with 
its fit from Figure [4] (magenta color). In the continuum limit, the mass of the 
non-Goldstone scPion will vanish and the fo state could become light close to 
the conformal window. The disconnected part of the correlator is required to 
resolve this issue. 

To study flavor-singlet mesons, we need to consider fermion 
loops which are disconnected (often called hairpin diagrams). 
Flavor-singlet correlators have fermion-line connected and 
fermion-line disconnected contributions from the hairpin dia- 
grams. To evaluate disconnected quark loops with zero mo- 
mentum, we need to sum over propagators from sources at each 
spatial location for a given time slice. To avoid the very costly 
0(V) inversions to compute all-to-all propagators in lattice ter- 
minology, random sources have to be used with noise reduction. 

A very interesting further challenge and complication is the 
existence of two types of distinct ++ scalar mesons. One of 
them is the composite fermion state and the other is the scalar 
glueball with the same quantum number. In dynamical sex- 
tet simulations, these two types of state will mix with an ob- 
servable spectrum of scalar mesons which will require a well- 
chosen variational operator set to disentangle the scalar state. 
This further underlines the room left for a light scalar state to 
emerge in the spectrum. It is also entirely possible that careful 
lattice calculations will shut down the Higgs interpretation. 

Staggered fermions present an additional complication from 
the contribution of pairs of pseudoscalar meson taste channels 
contributing to the scalar meson correlator. To be a physical 
state, the scalar meson fo has to be taste singlet. Taste selection 
rules then require that the fo meson couples only to pairs of 
pseudoscalar mesons of the same taste. We have shown earlier 
in Section 4 that the pion taste multiplet splits into the Gold- 
stone state and a variety of higher-lying non-Goldstone states, 
all degenerate with vanishing mass in the continuum limit. In 
the continuum limit only the taste singlet states (physical states) 
are expected to have the correct masses from the U(l) axial 
anomaly which is itself a taste singlet. The other non-singlet 
states remain light and create complicated threshold effects. 
This complication is present in the fo correlator masked by the 
physical two-pion intermediate state |94|. 



6.2. The Higgs particle and the dilaton 

If the sextet model is very close to the conformal window 
with a small but nonvanishing /3-function, a necessary condition 
is satisfied for spontaneous breaking of scale invariance gen- 
erating the light pseudo-Goldstone dilaton state. The model, 
as we argued earlier, is also consistent with chiral symmetry 
breaking (^SB) with the minimal Goldstone pion spectrum re- 
quired for electroweak symmetry breaking and the Higgs mech- 
anism. The very small beta function (walking) and^SB are not 
sufficient to guarantee a light dilaton state if scale symmetry 
breaking and ^SB are entangled in a complicated way. How- 
ever, a light Higgs-like scalar could emerge near the confor- 
mal window as a composite state, not necessarily with dilaton 
interpretation. To understand the important role of the non- 
perturbative gluon condensate in the partially conserved dilata- 
tion current (PCDC) relation and its related dilaton implica- 
tions, lattice simulations of the non-perturbative gluon conden- 
sate will be needed near the conformal window. 

For discussion of the PCDC relation constraining the proper- 
ties of the dilaton, we will closely follow the standard argument 
like in I69l l72ll73ll . We will also show how non-perturbative lat- 
tice methods can explore the implications of the PCDC relation 
when applied to the sextet model. 

In strongly interacting gauge theories, like the sextet model 
under consideration, a dilatation current & - &* v x v can be 
defined from the symmetric energy-momentum tensor &' v . Al- 
though the massless theory is scale invariant on the classical 
level, from the scale anomaly the dilatation current has a non- 
vanishing divergence, 

3^ = ©£ = ^G; v G^. (4) 

Although a(p) and G^ v G a/JV depend on the renormalization 
scale fi, the trace of the energy-momentum tensor is scheme in- 
dependent after renormalization. In the sextet model, the mass- 
less fermions are in the two-index symmetric representation of 
the SU(3) color gauge group. The gluon fields are in the adjoint 
representation with G" , a — 1,2, ...8. We will assume that the 
perturbative parts of the composite gauge operator G° v G fl/iv and 
0^ are removed in Eq. Q and only the non-perturbative (NP) 
infrared part will be considered in what follows. 

The dilaton coupling f a is defined by the matrix element 

im v {*)\<T{p)) = j(lf P v - g" v p 2 )e- ipx (5) 

with p 2 = m\ for the on-shell dilaton state cr{p). From the 
divergence of the dilatation current in Eq. Q we get 

(0\d M D"(x)\cr(p)) = Umle-v* . (6) 

The subtracted non-perturbative part of the energy-momentum 
tensor, 

[©£] = ^[ G "G £ H , (7) 

is defined by removing the perturbative part of the gluon con- 
densate in the vacuum, 

KL = P -J a 1 °% GmV - (OI^G^HOW . (8) 



The lattice implementation of the subtraction procedure will be 
briefly described after the derivation of the PCDC relation. 

It is easy to derive, like for example in [69], the dilaton ma- 
trix element of the energy-momentum tensor trace using some 
particular definition of the subtraction scheme, 

{cr{p = O)|[0^(O)] jvjj |O> - ^<0|[o£(0)] J0> . (9) 

When combined with Eq. (|6j, the partially conserved dilatation 
current (PCDC) relation is obtained, 



f- 

.1 (T 



1K(0)]J0> 



(10) 



Predictions for m a close to the conformal window depend on 
the behavior of f a and the gluon condensate G° v G fl/JV of Eq. 
There are two distinctly different expectations about the limit 
of the gluon condensate to f a ratio when the conformal win- 
dow is approached. In one interpretation, the right-hand side 
is predicted to approach zero in the limit, so that the dilaton 
mass would parametric ally vanish when the conformal limit is 
reached |69|. The formal parameter is the non-physical (frac- 
tional) critical number of fermions when the conformal phase is 
reached. In an alternate interpretation the right-hand side ratio 



of Eq. ( 10 1 remains finite in the limit and a residual dilaton mass 
is expected 117211731 . The two interpretations make different as- 
sumptions about the entanglement of ^SB and scale symmetry 
breaking but both scenarios expect a light dilaton mass in some 
exact non-perturbative realization of a viable BSM model. 

It is important to note that there is no guarantee, even with 
a very small /3-function near the conformal window, for the re- 
alization of a light enough dilaton to act as the new Higgs-like 
particle. Realistic BSM models have not been built with para- 
metric tuning close to the conformal window. For example, the 
sextet model is at some intrinsically determined position near 
the conformal window and only non-perturbative lattice calcu- 
lations can explore the physical properties of the scalar particle. 

6.3. The non-perturbative gluon condensate on the lattice 

Power divergences are severe in the calculation of the lattice 
gluon condensate, because the operator aG^ v G a,lv has quartic 
divergences. The gluon condensate is computed on the lattice 
from the expectation value of the plaque tte operator Up. On the 
tree level we have the relation 



lim fl _ 



1 



(1 -\ttUp) 



— { — GG)i att ice 

36 n 



(11) 



as the continuum limit is approached in the limit of vanishing 
bare lattice coupling go- At finite lattice coupling we have the 
sum of a perturbative series in go and the non-perturbative gluon 
condensate, 



a 6 ), 



(l - - tr Up) = V c„-gl"+a 4 - I l-GG) + 0( 

\ 3 H Zj *0 36 \P(g ))\n 'lattice v 

(12) 

where bo is the leading /^-function coefficient. There is no 
gauge-invariant operator of dimension 2 and therefore the or- 
der a 2 term is missing in Eq. ( 12 1. For small lattice spacing a, 



the perturbative series is much larger than the non-perturbative 
gluon condensate, and its determination requires the subtraction 
of the perturbative series from the high accuracy Monte Carlo 
data of the plaquette. The c„ expansion coefficents can be deter- 
mined to high order using stochastic perturbation theory [95 1. 
This procedure requires the investigation of Borel summation 
of the high order terms in the perturbative expansion since the 
coefficients c n are expected to diverge in factorial order and 
one has to deal with the well-known renormalon issues. The 
methodology has been extensively studied in pure Yang-Mills 
theory on the lattice 19611 . 

It will be very important to undertake similar investigations 
of the non-perturbative gluon condensate in the sextet model 
with full fermion dynamics. We hope to return to this problem 
in the near future. 

Summary and outlook 

We have shown that the chiral condensate and the mass spec- 
trum of the sextet model are consistent with chiral symmetry 
breaking in the limit of vanishing fermion mass. In contrast, 
sextet fermion mass deformations of spectral properties are not 
consistent with leading conformal scaling behavior near the 
critical surface of a conformal theory. Our new results are rec- 
onciled with recent findings of the sextet /3-function [3|, if the 
model is close to the conformal window with a very small non- 
vanishing yS-function. This leaves open the possibility of a light 
scalar state with quantum numbers of the Higgs impostor. The 
light Higgs-like state could emerge as the pseudo-Goldstone 
dilaton from spontaneous symmetry breaking of scale invari- 
ance. Even without association with the dilaton, the scalar 
Higgs-like state can be light if the sextet gauge model is very 
close to the conformal window. A new Higgs project of sex- 
tet lattice simulations was outlined to resolve these important 
questions. Plans include the determination of the S parameter 
and the sextet confining force with results on the string tension 
already reported, strongly favoring the^SB hypothesis [98|. 
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